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Alkali-earth atoms have a long-lived electronic excited state, and when atoms in this excited 
state are localized in the Fermi sea of ground state atoms by an external potential, they serve as 
magnetic impurities, due to the spin-exchange interaction between the excited and the ground state 
atoms. This can give rise to the Kondo effect. However, in order to achieve this effect in current 
atomic gas experiment, it requires the Kondo temperature to be increased to a sizable portion of 
the Fermi temperature. In this paper we calculate the confinement-induced resonance (CIR) for 
spin-exchanging interaction between the ground and the excited states of the alkali-earth atoms, 
and we propose that the spin-exchange interaction can be strongly enhanced by utilizing the CIR. 

We analyze this system by the renormalization group approach, and we show that nearby a CIR, 
the Kondo temperature can be significantly enhanced. 


I. INTRODUCTION 


Cold alkali-earth atomic gases have been widely used 
for building atomic clocks, with which the record of the 
most accurate optical lattice clock has been achieved [T]. 
This is because alkali-earth atoms have a very long-lived 
excited ^Pq state whose single-particle lifetime can be as 
long as many seconds. This excited state ^Pq and the 
ground state are viewed as two internal states of the 
“orbital” degree of freedom. Recently, there is an increas¬ 
ing experimental interest in studying many-body physics 
with alkali-earth atoms, including the SU{N) symmet¬ 
ric interaction and the orbital degree of freedom 0- 
[7]. In particular, recent experiments have demonstrated 
the inter-orbital spin-exchanging scattering between the 
ground state and this ^Pq state in fermionic ^^Sr [4] 
and ^^^Yb atoms EE]. 

Utilizing different AC polarizability of and ^Pq 
states, one can realize the situation that atoms in the 
^Po state experience a deep lattice and are localized, 
while atoms in the states experience a shallow lat¬ 
tice and remain itinerant, as shown in Fig. Therefore, 
due to the spin-exchanging scattering between these two 
states, atoms in the ^Pq state can play a role as magnetic 
impurities in the Fermi sea of atoms in the state, 
which can give rise to the Kondo effect [8] . Realizing the 
Kondo effect with cold atoms IBHIZ] can add a few new 
ingredients to the Kondo physics, such as the SU{N) 
Kondo model, manifestation of the Kondo effect other 
than transport properties and non-equilibrium dynam¬ 
ics. When the magnetic coupling is much weaker com¬ 
paring to the Fermi energy, such as in the cases of solid- 
state materials, the Kondo temperature is a few orders of 
magnitude lower than the Fermi temperature. However, 
with current cooling power, normally an atomic Fermi 


* Electronic address: pengzhang@ruc.edu.cn 
1 Electronic address: hzhai@tsinghua.edu.cn 





FIG. 1: Schematic of the system under consideration. The red 
balls are alkali-earth atoms in ^Pq state. They are trapped 
in a deep lattice and are localized, and have a lower density. 
The blue balls are alkali-earth atoms in ^So state. They are 
trapped in a shallow lattice and are itinerant. The arrows 
denote the nuclear spin degree of freedom. The system is 
confined into a one-dimensional tube. 


gas can only be cooled to ^ O.lTp (Tf denotes the Fermi 
temperature). Therefore, the most challenging question 
is how to increase the Kondo temperature to the range 
attainable by current experiments. 


In this paper we propose a scheme to overcome this 
challenge by using confinement-induced resonance (CIR). 
The CIR phenomenon describes resonant enhancement of 
the one-dimensional effective interaction strength when 
a system is confined into a quasi-one-dimensional tube 
nann. This phenomenon has been observed in previous 
cold atom experiments with alkali atoms [20]. Here we 
generalize the CIR phenomenon to the inter-orbital scat¬ 
tering between ^Sq and ^Pq states of alkali-earth atoms. 
We will show that a CIR can strongly enhance the spin¬ 
exchanging scattering, and consequently, the Kondo tem¬ 
perature can be increased to a sizable fractional of the 
Fermi temperature when a CIR is approached. 
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II. CONFINEMENT INDUCED RESONANCE 
A. Zero-Magnetic Field 

To illustrate the basic ideas, we first discuss a two-body 
problem at zero magnetic field. Let us briefly review 
interaction between two fermions in two different orbital 
states ^Pq (denoted by |e)) and ^Sq (denoted by |^)) and 
different nuclear spin states (for simplicity, here we only 
take two nuclear spin states denoted by | f) and | |)), 
respectively. We can introduce four antisymmetric bases 
for the internal degrees of freedom 

|±) = i(|5e)±|e5))(|n)TUt)) (1) 

|5t;et) = ^(|5e)-|e5))|n) (2) 

l5 4-;e4-) = ;^(|5e) - |e5))m)) (3) 

in which s-wave scattering is allowed. |+) is orbital 
triplet and nuclear spin singlet, and the other three are 
orbital singlet and nuclear spin triplet. Since nuclear 
spin does not participate the inter-atomic interaction 
process, the interaction part possesses the nuclear spin 
rotational symmetry, for which spin singlet and triplet 
will not mix and the interaction potentials are the same 
for |—) and \g t; e t) and \g e |) channels. Therefore, 
the inter-atomic potential U(r) is diagonal in the bases 

as 

V+{r)V+ + V-{y){V- + + V^): 

where Vi = |i)(z|, i = ±, and = \g \-,e \){g t;e t I 
and = \g e \) {g e |. The two interaction poten¬ 
tials are denoted by V±{r) = 27r/i^a±J(r)^(r-)//i {g is 
the two-body reduced mass), and a± are two independent 
scattering lengths. 

We can rotate the interaction potential i^(r) into an- 
other bases {l^r f,e i), |g 4,; e t), |5 t;e i)}, 

where l^r t; e J,) = (l/v^)(|+) + |-)) and \g 4-; e t) = 
(l/\/2)(|—) — 1+)), and V{r) becomes 

V: +V- V- -Vj. 

^ = - 2 - 2 - 

+ V.{Vn + Vu), ( 4 ) 

where = \g'\-,ei){g t;e4, |, V^^ = \gV,e^){g 4-;et |, 
and 5ex = |^ t5 ^ ^ t I- la the presence of a lat¬ 

tice as described in Fig. when atoms in the |e) state 
are localized as impurities while atoms in the \g) state 
remain itinerant, the off-diagonal (U_ — U+)/2 represents 
the process that an itinerant fermion exchanges its spins 
with impurities, and this spin-exchanging process is the 
essential process responsible for the Kondo effect [21]. 
Normally, this interaction strength is much smaller com¬ 
paring to the Fermi energy, and the Kondo temperature 
is exponentially suppressed m- 


Now let us consider atoms confined into a quasi-one- 
dimensional tube by a transverse harmonic trap. Here we 
consider the situation that the transverse confinement 
is the same for both and states, which can be 

achieved by applying a two-dimensional optical lattice 
in the x^-plane with the magic wave length [ 22 |. And 
we first consider the situation without lattice along the 
longitudinal ^-direction. The free-Hamiltonian Hq can 
be separated into the center-of-mass part and the relative 
motion part Pr, where 

= + + (5) 

For single interaction channel with a three-dimensional 
scattering length Ug, it is known that the interaction 
strength of the effective one-dimensional interaction po¬ 
tential go^{z) is given by [18] 

* = ( 6 ) 

g \ a±J 

w here ^ qq is the ground state wave function of Pr, a_L = 
^/gujjh is the harmonic length and C = 1.4603 • • • is a 
constant, g diverges when a± = which is known as 
CIR [18]. This resonance occurs when the energy of a 
bound state in the transverse excited modes matches the 
scattering threshold [19]. 

At zero-field, it is easy to show that |-h), |—), ^ t) 

and \g e 1) are all eigenstates of free Hamiltonian Hq. 
Hence, under confinement, these four scattering channels 
can be treated as independent channels, and the reduced 
one-dimensional interaction still takes a diagonal form as 

Rid = + g-'P- + goPn + goPu]^i^)^ C^) 

where g^ is related to a+, and g^^ go are related to 
a_ via Eq. respectively. Therefore, when ro- 

tated to {15 t;e i), [5 i; e t), 15 t; e t), Ifl e 4-)} 
bases, the spin-exchange interaction term is now given 
by {g- — g-^)6{z)/2. When a±_ Ca+, g^ diverges and 
g- remains finite; and when a± Ca-^ g- diverges and 
gj^ remains finite, as shown in Fig. [^a). Therefore, the 
spin-exchanging interaction becomes very strong nearby 
these two CIRs, and consequently, the Kondo tempera¬ 
ture can be dramatically enhanced. This is the basic idea 
of our proposal. 


B. Finite Magnetic Field 

In real experiment, there is always a finite magnetic 
held. Due to the difference in the Lande g-factor be¬ 
tween 1 ^) and |e) states, \g t; e 1) and \g e f) states 
differ by a finite energy 6 proportional to the magnetic 
held strength [23]. In another word, this leads to a mix¬ 
ing term between |+) and |—) given by J/ 2 (|-h)(—|+ h.c.) 
Els]. Therefore, these two channels can no longer be 
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FIG. 2: The one-dimensional interaction strength g- and 
gc(h^/(glO^ao) is taken as unit) as a function of a±/a~. (a) 
The magnetic field B = 0. Here pc = 0. (b) B = 35G. 

Here we use as an example and take S = 27r/ill2 x B, 

and therefore S = 25/iA;Hz at 35G. In both cases, go always 
behaves the same as p_ in case (a). For ^^^Yb we take a+ ^ 
10a- and a- = 200ao [24l [25]. The first CIR takes place 
around a± ~ Ca- and the second CIR takes place around 
a± ^ Ca+, though the exact locations of CIRs will be shifted 
by magnetic field. The insets in (b) show p’s nearby the first 
CIR and (5p = p_ — p+. 


treated independently, which brings out additional com¬ 
plication to our proposal. Under the condition 6 huj^ 
both two channels are retained in the one-dimensional 
effective mode. To deduce the one-dimensional interac¬ 
tional strength, our calculation follows the standard pro¬ 
cedure of CIR discussed before m, that is, one first ob¬ 
tains an effective one-dimensional scattering amplitude 
in both Ip 6 i) Id channels by solving three- 

dimensional Hamiltonian with the confinement potential, 
and then constructs an effective one-dimensional model 
which gives exactly the same scattering amplitude. 

The Hamiltonian for the relative motion of the two- 
body system with a confinement potential can be written 
as 


dimension, interaction does not mix |+) and |—) be¬ 
cause it respects the nuclear spin rotational symmetry. 
However, the finite Zeeman-held S in the single-particle 
Hamiltonian breaks the nuclear spin rotational symme¬ 
try. In the quasi-one-dimensional system, since the vir¬ 
tual processes to the transverse excited levels are taken 
into account when reducing dimensionality, the effect 
of this Zeeman energy term enters the effective one¬ 
dimensional interaction term Vib through the interme¬ 
diate state energy of these virtual processes. Hence, the 
effective ID Hamiltonian is expected to take the form as 

where 

ViB = [g^'P+ + d-'P- + dc{Sc + 5 ^) + do{Pn + 

( 10 ) 

It should be noticed that Vid no longer respects the nu¬ 
clear spin rotational symmetry, and contains a non-zero 
off-diagonal term pc that couples |+) and |—) channels. 
Similar effect has also been discussed in the CIR of spinor 
atoms [26]. p+, p_ and Pc depend on the strength of 
the Zeeman field, as shown in Fig. [^b). The detailed 
derivation of these parameters are given in the appendix. 
While po are still related to a_ via Eq. H since \g t; e t) 
and Ip l;e 1) states are still eigenstates despite of the 
presence of the Zeeman term. 

Secondly, when p+ diverges, p_ and pc will also diverge, 
as shown in Fig. |^b). Similar happens when p_ diverges. 
After the base rotation, the interaction term becomes 


= 


d+^d- 


+ dc ] + 


d++d- 


- dc Pit 


d- -d+ 


(Sex + + do (Pn + Pu) 


S(Z), (11) 


where the spin-exchanging term is still given by (p_ — 
p+)(5(z)/2. Fortunately, as shown in Fig. [^b), one can 
see that, nearby one of the CIRs, either p_ diverges much 
slower than p+, or p+ diverges much slower than p_. 
Therefore (p_ — p+)/2 still displays a divergent behavior 
and the insight gained from the zero-field limit will stay 
hold. 


III. THE KONDO EFFECT 
A. Lattice Model with a Single Impurity 


IIr = 


2 ^ . 2 
2 


(P++P_) + ^(5e + 


+ V+(r)P+ + V-(r)(P. + + Vu), 


Si) 

(8) 


where Sc = |+)(—|- Two important features at fi¬ 
nite S are worth emphasizing here. First, in the three- 


Now we consider turning on the lattice potential as 
shown in Fig. which localizes atoms in the |e) state as 
impurities in a Fermi sea of atoms in the |p) state. We 
also consider the regime where the density of impurity 
atoms is much more dilute than the density of itinerant 
atoms, and for simplicity, we consider a single impurity 
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problem. The tight-binding model is given by 

H = H^^H, (12) 

= cos k)cl^Cka + <5572 - T (13) 

/c,cr k 

= IE { 

/c,g ^ 

^Urikq^-^S^nkq-\-U2S%^ ^ (14) 


where c/^cr and da- are fermion operators for itinerant 
fermions and impurity fermion, respectively, 5'+ = 

S- = d\d^, = (l/2)(4c^t - ^kq = 4,iCg,t, 

4q = = (V2)(4,tCg,t - 4,iCgd) and = 

4 t^<?d ^ 4 ^ hopping amplitude of itiner¬ 

ant fermions. The interaction parameters J_l, Jz^ Uz, U-^ 
and U- are related to g- and gc via 

J± (X-{g^ - g-), Jz X-{g^ ^ g--2go), (15) 


U X 


+ ^- + 2^0 


Ui X -2gc, U 2 X gc, (16) 


where the proportional constant depends on the details 
of Wannier function overlap between localized atoms and 
itinerant atoms (In the plot of following figures, this con¬ 
stant has been chosen as 0.03). 

Here we first discuss a few features of this lattice model: 

1. The first three terms (J±-and J; 2 -terms) in Hi de¬ 
scribe the Kondo coupling. From Fig. one can see 
that i) in the left side of the first CIR, and the right 
side of the second CIR, J±,Jz < 0 and ii) in the right 
side of the first CIR, and in the left side of the second 
CIR, J±,J, > 0. Thus, both ferromagnetic and anti¬ 
ferromagnetic Kondo couplings can be accessed by tuning 
the confinement length a^. At zero-field, when go = g_^ 
from Eq. (15) we have J± = Jz and the Kondo cou¬ 


pling respects spin-rotational symmetry, while at finite 
Zeeman field, generally go is not equal to g- and it gives 
rise to an anisotropic Kondo model. 

2. The fourth term {U-teim) in Hi describes a poten¬ 
tial scattering that does not depend on spins. This comes 
from the diagonal term in Eq. (11). Eurthermore, at fi¬ 
nite field, due to the absence of the full spin rotational 
symmetry (a rotational symmetry along spin z-axes is 
still present), there exist other two scattering terms {Ui- 
and 1 / 2 -terms) in Hi. It naturally raises the question that 
whether these extra terms will affect the Kondo physics. 

3. The ^-term in Hq comes from the different ^-factor 
between localized and itinerant fermions. When this 
term is sufficiently large, it tends to polarize fermions 
and will destroy the Kondo physics. 

We shall also remark that here the interaction between 
itinerant fermions is ignored. Since microscopically it is 
described by another independent scattering length be¬ 
tween atoms in \g) states with different nuclear spins 
(normally denoted by Ugg.) The CIR for Ugg is reached 



FIG. 3: Second order diagram that can renormalize J-. Black 
arrows denote the itinerant atomic spin and red arrows denote 
the impurity atomic spin, (a) and (b) are “particle” process, 
(c) and (d) are “hole” process. 


at a different confinement radius when a± = Cugg, where 
the interaction between itinerant fermions will become 
very strong. With interactions, the effect of a magnetic 
impurity in a Luttinger liquid has been studied before 
imEH], and in the strongly interaction limit the results 
will be reported elsewhere m- Whereas, at these two 
CIRs we focus on here, this interaction is rather weak 
and can be safely ignore. 


B. Renormalization Group Studies. 


To address the effect of these extra terms and to make 
a more concrete predication of the Kondo temperature, 
here we adopt the renormalization group (RG) approach 
well established for the Kondo problem [21] |30] . The key 
idea is to iteratively integrate out the high-energy modes 
of itinerant fermions and see how the interaction param¬ 
eters flow. In this analysis, we only consider the lowest 
order virtual processes. For simplicity, only the scatter¬ 
ing processes that renormalize J_ are shown explicitly 
in Fig. as an example. The renormalization to the J_ 
term from the processes in Fig. [^a) when the cutoff is 
reduced from D to D — SD reads 




1 Jzl QZ t 

^ki-^p,i —^-W ^ 


i^S-S^Po\SD\ci^c,^^ 


- ^JziJ-Po\^D\S c\^Cqi 


(17) 


where po is the density of states of the itinerant atoms 
near Fermi surface and we have set CpcJ = 1 for the p 
states in the energy scale between D and D — 6D and 
= 1/2 for up-spin impurity. The approximation in 
Eq. ( [I7| ) comes from the fact that Cp ^ D and D ^ cj, e^. 
Similarly, one can write the renormalization contribution 
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to the J_ term from Fig. [^b),(c),(d) as 

-Jz^J-po\5D\S cl,^Cqi jj (IS) 

~ -^JziJ-Po\^D\S c\^Cq^ ]J ^ ^ (^^) 

-^Jz^J-Pq\5D\S c^^Cg4, —, (20) 

It turns out that the contribution of the diagrams which 
involve the 11^ term is zero. Thus one can sum up all the 
diagrams and obtain the renormalization equation for J_ 
as 


^ - j^t) 

dD 4 ^ 


D 


1 

~5^D 


( 21 ) 


By repeating similar procedures, one can find the renor¬ 
malization equations of all the interaction parameters in 
Eq. (14) as follows: 


dJ-\- pq 

—---j+j, 

dJ- _ j j 

dD 2 " 


- T T 

dD - 


dC 

dU2 _ Po J J 
-dD - 


1 1 

D-S D 
1 1 

D + S D 
1 1 


D-S 

1 


D + S 
1 


D-5 D + 5 


^=0 ^=0 
dD ' dD 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


where we have taken D as the energy cutoff, and its initial 
value is chosen to be an energy scale of the order of the 
Fermi energy. And for simplicity, we have also assumed 
a constant density of state denoted by po- When (5 = 0, 
Eq. (26) can reduce to the RG equations derived in Ref. 
[30] . At this level of approximation, U and Ui terms are 
renormalized. 

Here we focus on the antiferromagnetic case. In Fig. 
[^a), we show a typical flow of the RG equations with S 
Ep. We find that as lowering the energy cutoff, J± and Jz 
term diverge much faster than U 2 term. This in fact can 
already be seen in the RG equation, since dU 2 /dD scales 
with 6/D‘^ and therefore evolves very slowly when 6 
D. Thus, at the energy scale when J±, Jz diverge, the 
Kondo effect appears while the strength of other terms 
remain quite small. Hence, we conclude that the extra 
interaction terms in the lattice model Eq. (14) will not 
affect the Kondo effect in this system. 

The divergent energy scale for J± , Jz is normally taken 
as the Kondo temperature m- In Fig. gb) we show 
the dramatic increase of the Kondo temperature as ap¬ 
proaching one of the GIRs. The Kondo temperature can 
increase to ^ O.lTp which is attainable by current exper¬ 
iments. The underlying physics is basically the increas¬ 
ing of spin-exchanging coupling as we discussed above. 
In the plot of Fig. Sb) , we have restricted out initial 




FIG. 4: (a) Flow of parameters J±, Jz and U 2 as lowering en¬ 
ergy cutoff D. Here we have taken a typical density that gives 
rise to Ep = x 10^, and the initial values of J±, Jz and U 2 
as poJ± = 0.2377, poJz = 0.2376 and poUz = —3 x 10“^. (b) 
The Kondo temperature increases as the confinement length 
a± is tuned toward the first CIR from the anti-ferromagnetic 
coupling side. Here we have taken S = 0.07Ef. 


interaction parameter poJ±^poJz ^0-2. When it is very 
close to the CIR, the initial value of these interaction 
parameters are already very large which invalidates the 
perturbative RG approach. The impurity physics with 
a large or even divergent spin-exchanging interaction re¬ 
mains a challenging issue, and this makes the experimen¬ 
tal quantum simulation studies of this model even more 
interesting. 

Solving the RG equations for (5 ^ Ep or (5 Ep, we 
find none of the interaction parameters will diverge even 
when the cutoff D is lowered to zero as shown in Fig. 
which means the absence of the Kondo effect. Therefore, 
for a fixed a_\_/as, increasing 6 simply by increasing the 
magnetic field, it can drive a crossover from a Kondo 
regime to non-Kondo regime. 


IV. CONCLUSION. 

In summary, we have studied the CIR for two-orbital 
alkali-earth atoms with inter-orbital spin-exchanging in¬ 
teraction at finite magnetic field. We show that the CIR 
can strongly enhance the spin-exchanging scattering and 
hence dramatically increase the Kondo temperature. The 
signature of the Kondo effect will manifest not only in the 
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D/Ef 

FIG. 5: RG flow of parameters J±,Jz and in a strong 
magnetic field with B=100G and 6/Ef — 1.4 . 

transport properties, but also in other quantities such as 
spin susceptibility that can be measured by cold atom 
experiment [31]. Our proposal is very useful to ongo¬ 
ing experiments on many-body physics with alkali-earth 
atoms. 


Appendix A: Confinement-induced-resonance at 
Finite Magnetic Field 



e/s e/6 



FIG. 6: Scattering amplitude as function of incident energy. 
In (a) and (c): incident energy e ^ 0 and incident atoms are 
in state \g 1-; e t )(o = 0?/^ = !)• Here, a,/3 have been defined 
in Eqs. ( |A4|A5[ ). In (b) and (d): incident energy e ^ (5 and 
incident atoms are in state \g t; e l.)(a = 1, = 0). As shown 

in (a) and (b), 5 ~ 0.1/iCt;, the scattering amplitude given 
by ID model almost coincides with that given by complete 
calculation. While, if (5 ~ 0.4/ia;, the difference of scattering 
amplitude begin to show up, as shown in (c) and (d). 


To determin^+ , and gc^ we can re write the Hamil system the two-atom scattering wave function 

tonian in Eq. ^ as written as 


Hr = 


_2 1^^^ e 


V„ 


+ + Vi(>Sex + 


(^ti + 

(Al) 


|^(r)) = v|>(n )(^,^)|5 e + e (A3) 

where the functions (z, p) and are given 

by 


where it is not necessary to consider the \g t; e t) and 
\g e 1) channels at present since they do not cou¬ 
ple to other channels even in the presence of a mag¬ 
netic filed. Here the operators Vq and Vi are given 
by Vi = 2'K}T?asi6{v)^{r-)/ii {i = 1,2) with = 
(af -h <^7)/2 and a^i = {aj — af)/2. In Eq. (Al) we 
have shifted the threshold energy by a constant —6/2 . 
The relative motion of incident atoms are in t he trans - 
verse ground state (/)^=o,m;s=o(p): where p = , 

n = 0,l,2,-- - is the transverse principle quantum num¬ 
ber and rriz = n, n — 2, n — 4, • • • , is the quantum 
number for the angular momentum along the 2 )-direction. 
Here = 0 (1) when n is even (odd). Since the system is 
invariant under the rotation along the ^-direction, rriz is 
conserved and only the transverse states with rriz = 0 are 
involved in this problem. We further assume that both 
the Zeeman energy 6 and the relative kinetic energy e of 
the two atoms in the channel |p i; e t) are much smaller 
than the energy gap 2huo between the first transverse ex¬ 
cited state and the transverse ground state, i.e.. 


(A2) 


0n=O,m^=O (p) 

4^n,mz=o{p)’’! (''^4) 


^W>{z,p) = 

n=2,4,6,8,--- 

+ (A5) 


n=2,4,6,8,' 


Here we consider the general case where e could be either 
larger or smaller than 6. When e < 6, the incident atoms 
are in the state \g e f), which means P) = (0,1). 
When €> 6, the incident atoms are either in the state 
\g t; e |) or \g e f) which means either (a, /3) = (1,0) 
or {a^j3) = (0,1) for the system. In Eqs. ( |A4| A5 ) 
the function (j)n,mz{p) is the eigen-wave function of the 
transverse Hamiltonian (i.e. the Hamiltonian of a two- 
dimensional harmonic oscillator in the x-y plane with 
frequency uj) with quantum numbers {n^rUz), which sat- 
isfies (j)„^rn^=o{p = 0) = l/{^/na±) with a_L = ^h/{pu!) 
being the characteristic length of the transverse confine¬ 
ment. The parameters and are 


e, (5 ^ 2huj. 
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where the term huo in the right-hand-side of Eq. (A7) is 


FIG. 7: Red and blue solid line denotes the positions 
and of the first and second CIR, respectively. The blue 

dashed line denotes that if one extends the calculation to 
the limit of low trapping frequency, the position of the sec¬ 
ond CIR approaches the position of orbital Feshbach reso- 
nance(OFR) in 3D system [32], which is denoted by the red 
dashed line. The black dotted line is given by the condition 
a± = y/2h‘^/ii5. Our effective ID model is applicable only 
in the region with a± < jfi6, i.e. the region below the 


dotted line. 


given by 


^(tt) _ ^ 

l 2 fj,{e- 5 ) _ 

V " 


\l ’ 


' + - e) 




(A6) 


while the scattering amplitudes 

coefficients Bn'^^ and can be obtained from the 

Schrodinger equation 


contributed by the zero-point energy of the transverse 
ground state. 


Substituting Eqs. ( A4|A5 ) into Eq. (A7) and perform¬ 
ing the operation 


2 r-\-s POO 


on the both hands of Eq. (A7), we can obtain the rela¬ 
tions 


/(«)= 

. 20? 


(AS) 

R(tt) — 

-^n 

20r 


(A9) 

Kn'^'’aj_ 

fW) = 

. 20r 


(AlO) 

R(4.t) _ 

-^n 

20r 


(All) 

0t)„, ' 


where 


^(n) ^ ^ p = d) 

oz L 

^(tt) _ ^ z^f^^^\z, p = d) 
oz L 


2;^0+ 


z^0+ 


(A12) 

(A13) 


iIe|^'(r)) = (e + MI^'(r)), 


(A7) 


Eurth ermore, substituting Eqs. ( |A9|A11| ) into 
Eqs. ( A4|A5 ) and using the fact 0n,m7=oTp" = 0) = 
1 / (y^TTojjTwe obtain 






ae 




V^a_L 


I 3 e 


z 


V^a_L 






.2e 


|2;| 


A 




+ 


^ _ fk^ 
a± ^ I 2 




2e 




A 


ha±_ 


’ V 2 






ha±^ 


(A14) 


(A15) 


where the function A [^, z^] is defined as A [^, z^] = 

Vs' + u. As shown in Ref. [18] , this func¬ 


tion can be expanded as 


^ V — ^ + C 


2,1 + ^^ 


0(0, 


(A16) 


































































with being the Hurwitz-Z eta functio n. Substi- 

tutiM^qs. ( |A14|A15 ) into Eq. ( A12|A13 ) and using 
Eq. (A16), we can obtain the factors ^ and Eur- 

ther using these results and Eqs. (A8)-(A11), we can fi¬ 


nally obtain the expressions of the scattering amplitudes 
and 


-1 


i+iAP+ij:zA±^) 


' A'P 


(A17) 

where /, A, P and are e-independent 2x2 matrix 
with expressions 


I = 


1 0 
0 1 


0 


ry _ ( 

~ \ 0 k^A) 


(A18) 


and 


A = -^ 
2 


^sO ^sl 
^sl ^sO 



(A19) 


The expressions of A^ can be obtained by straightforward 
calculation we introduced above, and y e not needed in 
the following discussion. In Eq. (A17), ^ means the 
inverse matrix of [...]. 


On the other hand, according to Eqs. ( |A4 A5), in the 
long-range limit ^ oc (i.e., |^| ^ a±) the asymptotic 
wave function of the two-atom relative motion reads 

|^(r)) \g e |) 


+ 




0o,o(p)- 

(A20) 


Therefore, the long-range behavior of the two-atom 
relative wave function is completely determined by 

Eq. (IaTtI). 


Now let us consider a pure one-dimensional system 
with Hamiltonian Hid defined in Eq. Similar as 

above, in the bases {\g t; e i),\g i; e t)}, this Hamilto¬ 
nian can be re-written as 


+ 


9-\- P 9- P ^ . 9-\- P 9- ~ ‘^9c 

-^- Pn ^-^- 


+ 


9- -9+ 


(5ex + >Stj 


5{Z). 


^4,t 

(A21) 


With straightforward calculation, it is easy to find that 
the one-dimensional scattering wave function in this sys¬ 
tem is given by 


1^ 


ID 


{z)) = \g t; e |) 


+ 


\g g (A22) 


with a, /3, kPP defined the same as above, 

and the one-dimensional scattering amplitudes and 
are given by 


AU) 

Jm 

f(tt) 

^ID 


1 


/ + iAi^P 


(A23) 


where P is defined in Eq. (A18) and 
2h^ 


Aid = -- 


9 


9+ P 9- P 25^0 9- - 9+ 

9- -9+ 9+P 9-- 2^c 


(A24) 

is the one-dimensional scattering length matrix. 

By comparing Eqs. ( |A17|A19| ) with Eqs. ( |A23|A19| |, 
we find that in the systems where e is much smaller than 
hjuj so that the high-or der t erms i ( 2 ^ 0 ;) 

be neglected in Eq. (A17), when the one-dimensional 


parameters g± and gc satisfy 


f ^sO 

asi ^ 


V ^Sl 

j 



Ah 


^ (2 a + 2 ru ) . P . 

0 c(ii) 

5 + + 5- + 2(?c 9- — 9+ 

9- - 9+ 9++9-- 25 . 

(A25) 


we have the relation A zz Aid which gives 


fltt) 

'ID 


fin). fW) 


'ID 


fiW 


and 


|^(r)) <^„=o,„^.=o(p)|^lDW). 


(A26) 


(A27) 


Therefore, the Hamiltonian Hid in Eq. ([9|) with param¬ 
eters g± and gc given by Eq. (A25) is the correct one¬ 


dimensional model for these systems and the values of 
0 +, g- and gc at finite-magnetic field are plotted in Eig. 
I^b). In Eig. 1^ we p lot fiA) given by complete 

calculation (i.e. Eq. (A17)) a nd /| ]^ , from one¬ 

dimensional model (i.e. Eq. (A23)). It is shown that 
as long as S is small enough, the scattering amplitude 
given by one-dimensional model almost coincides with 
that given by complete calculation, no matter where the 
incident channel is located, as shown in Eig. [^a,b). How¬ 
ever when 6 is comparable to hw^ the difference between 
these scattering amplitude at high-energy regime begins 
to show up, as shown in Eig. [^c,d). 

In Eig. we show the show the characteristic length 
Ures^^ for the first and second CIR as function of Zee- 
man energy S. The position of the first CIR is insensitive 
to Zeeman field, while that of the second CIR is quit 
sensitive to Zeeman field. It is point ed t hat since our 
calculation is based on the condition ( |A2[ ), the effective 
one-dimensional model is applica ble only in the param¬ 
eter region S < 2huj or a± < ^J2hA jgS, i.e. the light 
region blow the dotted line in Eig. 
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